A FLARED, ORBITING, DUSTY DISK AROUND HD 233517 

M. Jura 

Department of Physics and Astronomy, University of California, Los Angeles CA 

90095-1562; jura@clotho.astro.ucla.edu 

O : ABSTRACT 

o 

(N 

We find that the infrared excess around HD 233517, a first ascent red giant, 
C/^ ■ can be naturally explained if the star possesses an orbiting, flared dusty disk. 

' We estimate that the outer radius of this disk is ~45 AU and that the total 

mass within the disk is ~0.01 Mq. We speculate that this disk is the result of 
the engulfment of a low mass companion star that occurred when HD 233517 
' became a red giant. 
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■ 1. INTRODUCTION 

2 ■ The evolution of dust disks into large sohd bodies such as asteroids, planets and comets 

eg ■ is of central importance in astronomy. While it is usually assumed that such disks are pri- 
mordial and associated with young stars, some binary systems such as the Red Rectangle 



apparently create orbiting reservoirs of gas and dust during their post main sequence evo- 
lution (see Waters et al. 1993, Jura & Kahane 1999). To date, there are only a few good 
candidate systems for orbiting dusty circumstellar material around evolved stars, and there- 
fore every example of this phenomenon is worth studying. Here, we propose that the puzzling 
infrared excess around HD 233517, a first ascent red giant, can be naturally understood as 
resulting from a flared, orbiting disk. 

Attention was originally drawn to HD 233517 because it is a K type star (my = 9.72 
mag) that is an IRAS source with an unusually large fractional infrared excess {Lm/L^ 
~ 0.06, Sylvester, Dunkin & Barlow 2001) produced by cold grains (T ~ 100 K), and it 
was thought to be an example of the Vega phenomenon - a main sequence star with dust 
(Walker & Wolstencroft 1988). However, detailed spectroscopic studies strongly suggest 
that HD 233517 is a first ascent red giant with a luminosity near 100 Lq (Fekel et al. 1996, 
Balachandran et al. 2000, Zuckerman 2001). While infrared excesses are common around 
second ascent red giants on the Asymptotic Giant Branch with luminosities in excess of 
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1000 Lq (see, for example, Habing 1996), relatively few first ascent red giant stars with their 
lower luminosities have detectable infrared excesses (Judge, Jordan & Rowan- Robinson 1987, 
Zuckerman, Kim & Liu 1995). HD 233517 has a uniquely high value of L/^j/L* among first 
ascent red giants where this ratio is generally less than 10^'^ (sec Zuckerman et al. 1995). It 
is unlikely that HD 233517 is a pre- main sequence star since it does not lie near any known 
region of star formation and since it lies far from the location of young stars in the H-R 
diagram (Fekel et al. 1996, Balachandran et al. 2000). Also, HD 233517 appears to be 
similar to a few other K giants which rotate rapidly, have strong lithium lines, a detectable 
far-infrared excess and are definitely post main sequence since at least in the case of PDS 
365, the carbon isotope ratio, ^^C/^^C, is approximately 12 (Drake et al. 2002). 

Three models to explain infrared excesses around first ascent red giants are (1) dust is 
being produced by mass loss, (2) nearby interstellar dust is illuminated accidentally by the 
red giant, and (3) the dust is orbiting (see Jura 1999, Kalas et al. 2002). The fraction of 
red giants that display true dust excesses is uncertain and controversial (Jura 1990, Plets et 
al. 1997, Jura 1999, Kim, Zuckerman & Silverstone 2001). In any case, since ground-based 
imaging shows that the 10 /im excess is physically associated with HD 233517 (Skinner et 
al. 1995, Fisher et al. 2000), at least for this particular star, the dust producing the infrared 
excess either is orbiting or is expanding away from the star. For some other stars, the infrared 
excess has been resolved with ISO and in these objects, the dust is probably not orbiting 
(Kim et al. 2001). 

There arc difficulties with the model that HD 233517 is currently losing enough mass 
to produce the observed infrared excess. As with other first ascent red giants with dust, 
the infrared spectrum of HD 233517 peaks near 60 //m. This spectral energy distribution is 
characteristic of cool material and is very different from that associated with a continuous 
mass loss rate where there is a substantial amount of dust near a temperature of 1000 K (see, 
for example, Sopka et al. 1985). One possibility is that the mass loss rate from HD 233517 
is episodic and currently the star is not expelling much matter. However, as noted by Jura 
(1999), the characteristic time scale for the dust to expand around a first ascent red giant 
to its inferred location may be less than 20 years. There is no evidence for recent infrared 
variability of this star, although the star does exhibit a full amplitude optical variation of 
0.02 mag with a 47.9 day period which may be caused by star spots rotating in and out of 
the view (Balachandran et al. 2000). 

Since the dust around HD 233517 may not be carried in a wind, we consider models 
where the dust is orbiting the star. One scenario is that the dust around HD 233517 is simply 
"left over" from the main sequence phase. However, a major difficulty with this model is 
that the fractional infrared excess around HD 233517, L/r/L*, of 0.06 is much larger than 
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that found for even the "dustiest" main sequence stars such as (3 Pic and HR 4796 where 
Lir/L^ is 2-4 X 10~^ (see, for example, Zuckerman 2001). 

Here, we speculate that HD 233517 had a low mass companion which was engulfed when 
it became a red giant. Although uncertain, the engulfment of this companion could have led 
to the ejection of an equatorial ring of orbiting material (see Taam & Sandquist 2000, Spruit 
& Taam 2001). As discussed by Pringle (1991), after this ring is created, the matter expands 
under the action of tidal torques which transfer angular momentum from the binary into 
the circumstellar material which thus evolves into an extended disk. Although not required 
in Pringlc's model, it is possible that in this dense equatorial ring, dust grains formed. We 
suggest that the structure of the circumstellar system is described by the models of Chiang 
& Goldreich (1997) for passive, orbiting disks which flare in response to the illuminating 
radiation. Below, we present details of this idealized model. 



2. THE DISK AND THE STAR 

2.1. Stellar parameters 

First, we list our assumed parameters for the star. Following Balachandran et al. (2000), 
we adopt a distance, of 620 pc, an effective temperature, T*, of 4475 K a radius, i?*, 
of 1.1 X 10^^ cm and a luminosity of 100 L0. We assume that the mass of the star, M*, is 
approximately 1 Mq, consistent with the star's location on the H-R diagram and measured 
surface gravity (Balachandran et al. 2000). The rotational velocity of the star is vsini = 
17.6 km s~^ which implies if the star rotates as a solid body and if the moment of inertia 
of the star is 0.1 M* Rl (taken from Schwarzschild 1958) that the angular momentum of the 
star is at least 4 x 10™ g cm^ s~^. 



2.2. Predicted Infrared Emission 

Dust is required to explain the emission from HD 233517 for A > 10 fim (see Skinner 
et al. 1995). Here, we use the model of a opaque, flared disk which is locally in vertical 
hydrostatic equilibrium (Chiang & Goldreich 1997). In this model, the temperature, Tdisk, 
as a function of distance from the star, R, is (see Jura et al. 2002): 
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where ks is Boltzmann's constant, G is the gravitational constant and fi is the mean molec- 
ular weight of the gas. Here, we assume that the gas is primarily H2 with [He]/[H2] = 0.2 
so that n = 3.9 x 10~^^ g. In view of the many uncertainties, and to be consistent with 
inferences for disks around pre-main sequence stars (Chiang et al. 2001), we assume a simple 
model where the half-thickness of the vertically-isothermal disk is taken equal to the param- 
eter, h, in Jura et al. (2002) which is a factor of \/2 larger than h in Chiang & Goldreich 
(1997), that characterizes the Gaussian density distribution above the plane. 

Since the disk is opaque, we can determine the flux from the source, F^,, by: 



= B,{Ta.sk)RdR (2) 



where Rout denotes the outer boundary of the disk where the temperature is Tout- For 
simplicity, we assume a sharp outer boundary although an exponential decay would be more 
realistic (see Jura et al. 2002). It is convenient to introduce the dimensionless parameter, x: 

From equations (1) and (2): 

287r cos i Rl [ keT, \ {kB%f [ 2 keT^R, \ ^'^ T^* x^^'^ 



Dl \ hu J {hcf \AQGM^n 



\ ^/•^ rxout „ll/3 



At high frequencies with Xout > 10, the observed flux is insensitive to the outer boundary 
condition since the integral in equation (4) is approximately 15. The spectrum is predicted 
to vary as u~^/^. At low frequencies where Xout <2, we may re- write equation (4) to find 
that: 

28n cosiRl f \ ^^'^ ( 2 fc^T^i?, ^ '^'^ 



If the disk is opaque, then at low frequencies, Fj, varies as u^. 

For the data in the IRAS Faint Source Catalog, between 12 //m and 60 //m, F,^ varies as 
in good agreement with the prediction that should vary as Furthermore, 
with Xout >> 1 and cos i = 0°, from equation (4), the predicted value of Fj^(25 fj,m) is 
3.9 Jy which agrees with the observed value of 3.4 Jy. Thus, without any "tweaking", the 
predictions are close enough to the data that the model should be taken seriously. 

Near 60 /im, the slope of the spectral energy distribution turns over. In the flared disk 
model, the peak in the spectral energy distribution is governed by the outer boundary of the 
disk. As shown in Figure 1, we have found that a model with Tout — 70 K, which occurs at 
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Rout — 45 AU, can reproduce much of the infrared data. However, at 1.35 mm, much more 
flux is predicted than the observed value of 2.6 ± 1.0, which is conservatively interpreted as a 
3(7 upper limit of 3 mJy by Sylvester et al. (2001). A possible resolution of this discrepancy 
is discussed below. 



2.3. Predicted Millimeter Continuum Emission 

The flared disk model drastically overestimates the flux at 1.35 mm. It is observed that 
between 1.35 mm and 100 /xm, Fj, varies as z/^'^ which suggests that at least at the longer 
wavelength, the system is optically thin. If the circumstellar dust particles are smaller than 
100 in diameter, then the opacity at 1.35 mm is probably considerably smaller than the 
opacity at 100 fim. Although the composition and size distribution of the dust particles arc 
unknown; we adopt parameters representative of circumstellar "silicates" from Ossenkopf, 
Henning & Mathis (1992) and Pollack et al. (1994). Thus, we take x^(100 jim) = 25 cm^ 
g~^, and for A > 100 //m, we assume that Xu scales as u'^ so that x,^(1.35 mm) = 0.14 cm^ 

To calculate the flux for the optically thin regime, we need to describe the mass dis- 
tribution. If E denotes the total dust surface density of the disk, then we approximate the 
model of Pringle (1991) for a circumbinary disk by assuming that all the matter is confined 
within angular radius (pout, where (f) — R/D^, and that within this region, 

E = Eo {^^y (6) 
With this description, if M^ust denotes the mass of dust in the disk, then: 

~ 47rR2 

^ -'^out 

The flux from the disk, F^isk, is: 

/(pout 
2 7: (t>d(t>B,{T) (1 - exp-^^^'^-'^^'^^/'^"'') (8) 


We evaluate equation (8) for the condition that at i? = Rout (or 45 AU), T^out = 0.04 g cm~^. 
As discussed below, this value of Equ* is chosen so that the disk is opaque at 100 /xm at its 
outermost region. The results are shown in Figure 1. With the revised model, we see that 
the predicted flux at 1.35 mm is consistent with the upper limit reported by Sylvester et al 
(2001). 
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2.4. Predicted CO Emission 

The flared disk model requires that gas be in hydrostatic equihbrium vertical to the 
orbiting plane of the system. We can estimate the expected amount of CO emission from 
HD 233517 if the gas and grain temperatures are equal to each other. First, 

T^T^tl^^y (9) 

If the gas is optically thick, then in a telescope of beam size, fltei, the observed brightness 
temperature, T^bs is 

-^"^ 27r0T(0) # ^ 14 7rT,„,0Lcog^ ^^^^ 

^tel 11 ^tel 

Within the IC/.'S beam of the IRAM telescope. Jura & Kahane (1999) placed a la upper limit 
to the CO J 2— >1 emission line of 29 mK. In our model, (pout = 0'.'073, and the predicted 
brightness temperature in this line is 17 mK. Thus the model is consistent with the observed 
upper limit to the CO line emission. 




2.5. Disk Mass and Angulcir Momentum 

We derive a minimum mass of the disk from the criterion that our model requires the 
disk to be opaque even at it! = Rout at 100 /xm or T,out Xj/(100 /im) > 1. With Xi,(100 iim) = 
25 cm^ g~^, this imphes that T,out > 0.04 cm^ g~^. Therefore, from equation (7), the mass 
of the dust in the disk is 2 x 10^^ g. If the gas to dust ratio is 100, the total disk mass is 
~0.01 M0 

Another important parameter is the disk angular momentum, J. With a sharp outer 
boundary and the surface density variation given by equation (6), we may write: 

J = ^^GM.Rout (11) 

We therefore estimate a total angular momentum of the disk of 3 x 10^^ g cm^ s~^, which 
is larger than the total angular momentum of the star estimated above to be > 4 x 10^° g 

2 —1 

cm s . 



3. DISCUSSION 



We have found that the infrared fluxes from HD 233517 are naturally explained with the 
simple flared disk model. The origin of this disk is uncertain. One possible model to account 
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both for the creation of a disk and the unusually high rotation speed of the red giant is that 
while on the main sequence, HD 233517 had a companion at an orbital separation of about 
1 X 10^^ cm (see Soker, Livio & Harpaz 1984). We speculate that when HD 233517 became 
a red giant, it engulfed this companion with the consequences of both creating an orbiting 
ring of circumbinary matter and spinning up the primary star. With the parameters derived 
above, the total angular momentum of HD 233517 and its associated disk is 3 x 10^^ g cm^ 
s~^. This amount of angular momentum could have been supplied by a companion star of 
0.12 M0 in a circular orbit of radius 1.1 x 10^^ cm. 

HD 233517 is a lithium-rich star (Fckcl ct al. 1996, Balachandran ct al. 2000). While 
single-star models might explain this result (de la Reza et al. 1997, Jasnicwicz ct al. 1999, 
Charbonnel & Balachandran 2000), in our picture, this lithium enrichment results because 
the system was a binary. The absence of ^Li argues against the simple idea that incorporation 
of lithium from a planet or low mass companion into the envelope of the star (Balachandran 
et al. 2000) can explain the observed abundance of this clement, but the capture of a 
companion might have led to an episode that stimulated interior mixing of lithium to the 
surface (Denissenkov & Weiss 2000). 

With a value of L/r/L* > 10~^, HD 233517 is a distinctive red giant. In the sample of 
Zuckerman et al. (1995) of more than 40,000 of the nearest luminosity class HI giants, the 
highest value of Lir/L^ is 10^^. Although the statistics are poor, perhaps only 10"*^ of all 
first ascent red giants have as much circumstellar opacity as docs HD 233517. We suggest 
that HD 233517 is unusual for two reasons. First, according to Duqcnnoy & Mayor (1991), 
only ~2% of all solar-mass main sequence stars have detectable companions with orbital 
periods near 10 days. Such binaries are the kinds of systems that can evolve into a red giant 
with a substantial disk as is found around HD 233517. Furthermore, once the disk is formed, 
it must have a limited lifetime. The massive dusty disk around HD 233517 may evolve in a 
fashion somewhat like disks around pre-main sequence stars. For A type stars which have 
luminosities somwhat less than that of HD 233517, the lifetime of the disk as a system with 
Lir/L^ > 10~^ is near 3 x 10® years, although some small amount of dust could linger for 
as long as 4 X 10^ years (see Habing et al. 1999, Spanglcr et al. 2001). Since the lifetime of 
a 1 M0 red giant as it evolves from 10 Lq to its maximum luminosity of 2700 Lq is ~2 x 10^ 
years (Girardi et al. 2000), only ~ 1% of the disks around red giants might have Ljr/L^ as 
large as 10~^. Given these two considerations, it is not too surprising that a system hke HD 
233517 occurs among first ascent red giants with a frequency of perhaps 10~^. A plausible 
explanation for the inferred decrease of with time is that the dust particles coalesce 

into larger bodies. 
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4. CONCLUSIONS 

We have modeled the currently available infrared data for HD 233517 and find that the 
observations can be naturally explained by a fiared, orbiting disk of mass ~0.01 M© and 
outer radius of ~45 AU. We speculate that this disk was created by the engulfment of a low 
mass stellar companion. 

This work has been partly supported by NASA. 
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Fig. 1. — A plot of the model spectral energy distribution for HD 233517 compared to 
observations. The squares show the data from IRAS; the upper limit to the flux at 1.35 mm 
(Sylvester et al. 2001) is shown as the upside triangle. The solid line shows the model for 
the opaque, flared disk; the dashed line shows the modiflcation to this model if the outer 
portion of the disk becomes optically thin at A > 100 yum as described in the text. 



